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We present a self-consistent mean field model of the extraction of atoms from a Bose-Einstein 
condensate to form a CW atom laser. The model is based upon the Hartree-Fock Bogoliubov 
equations within the Popov approximation, modified by the inclusion of spatially dependent source 
and sink terms, which lead to current flow within the condensate. The effects of this current flow are 
investigated for traps containing Rubidium (repulsive effective interaction) and Lithium (attractive 

qq ■ interaction) atoms. The extra kinetic energy associated with this flow is shown to be able to stabilise 

Q\ \ the condensate in the attractive case against mechanical instability. 
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The possibility of producing an intense coherent beam of atoms from a trap containing a sample of Bose-Einstein 
condensed atoms, a so called atom laser, has attracted much interest in the time since Bose-Einstein condensation 
was first achieved experimentally in a dilute gas of atoms jjj. To our knowledge little or no work has appeared 
however, investigating what effects extracting such a beam of atoms would have on the properties of the condensate. In 
particular, does the extraction of a small proportion of the condensate over some reasonable timescale drastically alter 
properties such as the condensate profile and noncondensate fraction 0? Given the great qualitative and quantitative 
success of mean field theories |J , specifically Hartree-Fock Bogoliubov (HFB) theory jij , in determining the collective 
excitation spectrum, condensate profile, noncondensate density and noncondensate fraction for trapped Bose gases, 
^v^j ' we develop a version of the HFB formalism within the Popov approximation j^] (which neglects two-body correlations) 
y—{ \ P| , in which we include a spatially dependent sink term, which extracts atoms from the condensate. This corresponds 
to the output coupling of atoms from the trap in some manner. We envisage some form of continuous Raman output 
coupling Q], but the specific mechanism is unimportant in this treatment. Also included is a stimulated pumping 
term which corresponds to the Bose enhanced scattering of atoms from the non-condensate into the condensate. The 
total number of atoms in the trap is kept constant and the noncondensate and condensate profiles are determined 
self-consistently. This corresponds to an assumption that the trapped atoms are replenished at the same rate as atoms 
are removed from the condensate via a thermal bath of atoms at a given temperature and that these 'new' thermal 
atoms equilibrate on a timescale that is small compared to the characteristic timescale of the trap. 
■ This leads to a modified time-dependent Gross-Pitaevskii equation now given by, 

e . 

m ^M. = j-|^V 2 + V ext (v) + g[n c (r) + 2n(r)] - | 7c (r) + f r«(r)} *(r, t). (1) 



Here we have made the usual decomposition of the Bose field operator, ip(r,t), into condensate and noncondensate 
parts, i.e., V>(r, i) = ^(r, t) + ■ip(r, t). The terms involving the interaction strength, g — AirTi 2 a/M, arise from the 
use of a contact interaction, g8(r), where a is the scattering length measured for binary scattering in vacuo. The 
condensate and noncondensate densities are given by n c and h respectively. The factor of 2 before the n(r) comes from 
the equivalence of the direct and exchange energies in the Bose case when one uses a contact interaction. The term 
involving 7 c (r) is the sink term due to the extraction of condensate atoms via the output coupling into the atom laser 
beam, with a coupling strength determined by j c . The Bose enhanced pumping of atoms from the noncondensate 
into the condensate is given by the term involving T, which is a constant determined self-consistently such that the 
number of atoms scattered from noncondensate to condensate per unit time is equal to the number of atoms extracted 
from the condensate. Note that this form of the coupling automatically includes the Bose enhancement factor. 

We now look for a steady state solution for the self-consistent condensate and non-condensate profiles by looking for 
solutions of the form ^(r, t) = <f>(r)e _l K'e l6l ( r ). Here /i is the chemical potential and 6(r) is a position dependent phase. 
Normally the phase is an arbitrary constant and is set to zero without loss of generality. The phase only becomes a 
relevant quantity in the context of phase differences or gradients. In this problem, with spatially dependent sink and 
source terms for condensate atoms, one does have a phase gradient which corresponds to a flow of atoms within the 
condensate. 

Making this substitution in the time-dependent Gross-Pitaevskii equation and separating into real and imaginary 
parts, one obtains two equations. The first is a modified time-independent Gross-Pitaevskii equation; 
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g[n c (r) + 2n(r)]}$(r)= fi^r), (2) 



the second being the differential equation determining the phase gradient 

#w +2 w*i + 2 *w «w + 1 (r - (r) _ 7c(r)) _ (3) 

Here we have made the additional simplifying assumption that the trap is isotropic. 

Similar analysis leads to the usual set of coupled HFB equations (within the Popov approximation) given by 

Cui{r) - gn c {r)vi{r) = EiUi(r) (4) 
Cvi{r) - gn c (r)ui(r) = -E t v t (r), (5) 

where the Hcrmitian operator, C, is modified by the addition of a kinetic energy term associated with the flow of 
condensate atoms and is given by 

L = ~&i v2 + &i irW) 2 + Vext{r) + 2g{nc{r) + " (r)) ~ (6) 

The HFB equations define the quasiparticle excitation energies and amplitudes Ui and Uj. The condensate density 
is defined such that, n c = |<3?(r)| 2 and noncondensate density such that, 

fi(r) = (V^(r)^(r)) = £ {k(r)| 2 + [\ Ul (r)\ 2 + h(r)| 2 ] N(E t )} , (7) 

i 

where N(Ei) is the usual Bose factor with a temperature determined by that of the external bath of atoms. 

These now form a closed set of equations which can be solved self-consistently giving the condensate and noncon- 
densate densities, excitation spectrum and phase gradient. 

We now present results for two different regimes. Firstly we consider an isotropic trap with V ex t(r) = V ex t(r) — 
±Afw 2 r 2 , where M( 87 Rb) = 1.44 x 1CT 25 kg, and the trap frequency v = lu /2tt = 200 Hz. The scattering length 
for &7 Rb is positive (repulsive effective interatomic interaction) and is given by a ~ llOao = 5.82 x 10~ 9 m. We 
consider a trap containing 4000 atoms in the steady state at a temperature of 90 nK. Solving self-consistently for 
the noncondensate fraction with no output coupling yields 2263 atoms in the condensate and 1737 atoms in the 
noncondensate at this temperature, i.e. just less than 50 % thermal depletion. We then extract condensate atoms in 
a localised region at the centre of the trap via a coupling j c (r) = Ae~ ar . The total number of atoms extracted in a 
characteristic time r = 1 /loq is given by, 



N nut = 



J d\n c {r)[\-e-^ T ]. (8) 



We choose a — 2 such that the full- width half-maximum (FWHM) of the coupling is of the order of 1/3 of the 
FWHM of the condensate. This enables one to induce moderately large phase gradients (and hence current densities, 
J = n c d9/dr) within the condensate. If one is considering a Raman output coupling scheme, the FWHM would 
correspond to the focusing of the laser beams. This is of course limited ultimately by the wavelength of the light used. 
With the parameters used here the focusing would have to be on the order of 1 /xm, which is on the extreme limit 
of achievability. We would like to emphasise that we are looking at a trap containing only 4000 atoms (to which we 
are limited numerically). Realistic traps are much larger now || and hence the focusing in experiments can be much 
weaker and attain similar results qualitatively to those presented here. 

In Fig 1 we present results for the condensate and noncondensate densities and the current density for a range of 
amplitudes, A, of the coupling strength. Note here that the current density is defined such that current flow towards 
the centre of the trap is positive. For relatively weak coupling the phase gradient established is small and the density 
profiles of the condensate and noncondensate are essentially identical to those obtained without any laser output. With 
increasing coupling strength, the current densities grow monotonically (left panels) with a corresponding increase in 
the effects upon the densities ||. The condensate becomes slightly depleted at the centre and the noncondensate 
density increases slightly. This can be thought of as an effective heating caused by the extraction of the laser beam, 
but even in this extreme case where we have tried to maximise the phase gradients by taking a small extraction region, 
the total increase in the number of atoms in the noncondensate as opposed to the condensate is from 1737 atoms with 
no extraction to 1822 atoms at the maximum extraction rate when N out ~ 100 atoms in time r. 
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When the coupling strength is increased further (right panels) the effects upon the densities become more pro- 
nounced and the condensate density at the centre of the trap becomes depleted significantly. This actually leads to a 
reduction in the rate of output of atoms from the trap as the coupling strength is increased. This can also be seen 
from the fact that the current density for the case with strongest coupling (dashed lines), has actually decreased. 
There is therefore an optimal strength with which to couple atoms from the condensate so that one can maintain a 
reasonable yield of atoms without depleting the condensate density excessively. This is equivalent to saying that there 
is a maximum flow rate for atoms within the condensate. 

The second case we wish to consider is the case when the trapped atoms have a negative scattering length, i.e. 
attractive effective interactions. We therefore consider a trap with i/q = luq/2ti = 150 Hz and M( 7 Li) = 1.16 x 10~ 26 
kg. The s-wave scattering length for 7 Li is given by a = —1.44 x 10 -9 m. The parameters here correspond closely to 
those of the Rice group jl0| . 

It is well known that condensates with attractive interactions are only metastable at best and above some critical 
condensate number become unstable to mechanical collapse, which is characterised by the I = collective mode going 
soft, i.e. the lowest / = excitation frequency goes to zero (0]]. The critical number for the above parameters has 
been shown to be N c = 1241 atoms at T = 0. It has also been shown that increasing the temperature decreases 
the stability and decreases the critical number N c . This decrease manifests itself in a decrease of the / = mode 
frequency with temperature. For example, with 1150 atoms in the condensate and a total number of atoms in the 
trap determined self-consistently from the HFB equations, the I = mode frequency drops from u>i = o — 1.523u;o at 
T = to cl>; =0 = 1- 428wo at T — 75 nK, with 1312 atoms in the noncondensate, viz. a total of 2462 atoms in the trap. 

The (meta) stability of the condensate is due to the kinetic energy associated with the trapping potential - there 
is no stable condensate in the uniform case with attractive interactions. The extraction of atoms via an atom laser 
sets up a current flow within the condensate which has the effect of increasing the kinetic energy. It is therefore 
feasible that the extraction of atoms via some method of localised output coupling could stabilise the condensate 
against mechanical collapse. We therefore apply our model to the case with 1150 atoms in the condensate at 75 nK. 
In this case a = 256 which corresponds to a FWHM of 0.4 //m. This puts one into the uv region for Raman output 
coupling, which is infeasible. However for this trap a very small output region is required as the characteristic size 
of the condensate is small. To experimentally investigate these effects one would need to use a trap with a smaller 
confining potential and hence larger (less dense) condensate. The amplitude of the coupling is 2500 ujq which yields a 
cw laser output of 129.46 atoms/r. The condensate, noncondensate (magnified by a factor of 20) and current densities 
are shown in Fig 2. The peak in the current density approximately corresponds, in terms of position to that of that 
of the point of inflection in 7 c (r). The relative energy scales of the kinetic energy due to the current flow and the 
scale of the interaction energy in this case is of the order 



increased form u)i—q = 1.428wo to ll>i = q = 1.464cjo- This represents a significant decrease in the softening of the mode 
and a strong indication of the feasibility of stabilising the condensate against collapse. 

In conclusion, we have developed a model based upon well tried self-consistent mean field theory that incorporates 
the effects of continuous extraction of atoms from the condensate and pumping via scattering of atoms from the 
noncondensate to condensate. We have investigated the effects of the induced current within the condensate upon 
the steady state properties of both the condensate and noncondensate. For trapped atoms with repulsive interactions 
and moderate extraction rates, the effects were not found to be excessively damaging and it appears likely that atoms 
could be extracted by means of some form of output coupler in a continuous manner, without significantly effecting 
the steady state properties of the condensate. 

In the case of trapped atoms with attractive interactions the spatially localised extraction of an output beam of 
condensate atoms, setting up a flow of current within the condensate was found to increase the kinetic energy of the 
condensate. This was found to reduce the softening of the / = mode through thermal effects, hence stabilising the 
condensate against mechanical collapse. It is hoped that for larger traps (weaker confinement) larger phase gradients 
may be created which, further increasing the kinetic energy, may be able to increase the critical number of condensate 
atoms for given trap parameters beyond that of the T — case with no extraction. 
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acknowledge the Royal Society (London), University of Auckland Research Committee, the Marsden Fund of the 
Royal Society of New Zealand and the UK EPSRC for funding. Thanks go to Tony Wong, Dan Walls and Keith 
Burnett for many useful discussions and very special thanks in general to the quantum optics group at the University 
of Auckland for their warm hospitality. 
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Even with this relatively small addition to the kinetic energy the frequency of the I = collective excitation has been 
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FIG. 1. Condensate and noncondensate densities in units of d 3 and current densities (ujod 2 ) for a trapped gas of 4000 
87 Kb atoms at T = 90 nK with a spatially dependent output coupling rate given by 7 C = Ae~ ar . a = 2. Left panels; A = 10 
(solid line), A = 80 (dotted line), A = 120 (dashed line). Right panels; A = 160 (solid line), A = 200 (dotted line), A = 240 
(dashed line). 



FIG. 2. Condensate, noncondensate and current densities for a trapped gas of 7 Li atoms at T = 75 nK with 1150 atoms in 
the condensate and Gaussian output coupling strength 7 c (r) centred at the origin. 
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